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We show that a "fragmented Bose condensate" in 
which two or more distinct single-particle states 
are macroscopically occupied by the same species 
of boson is inherently unstable to the formation 
of a conventional Bose condensate whose macro- 
scopically occupied state is a linear combination 
of the "fragments" with definite relative phases. 
A related analysis shows that a reproducible rela- 
tive phase develops when two initially decoupled 
condensates are placed in contact. 

The essential property of a Bose condensate is the ex- 
istence of a single macroseo^cally large eigenvalue of its 
one-body density matrix. EJ~I3 This property reveals the 
existence of long-range phase coherence across the fluid, 
which leads to the two-fluid description of Bose conden- 
sates and implies the quantization of circulation in su- 
perfluid flow.aB Hypothetical Bose states without this 
property were discussed long ago by Nozieres and St. 
JamesErQ, and have recently received renewed theoretical 
attention in the case of trapped bosons.BB In particular, 
Wilkin et alB found that the azimuthally symmetric ro- 
tating states of a dilute Bose gas with attractive interac- 
tions are "fragmented" Bose condensates in which many 
single-particle states (of differing angular momenta) are 
macroscopically occupied. Such multi-condensate Bose 
fluids would be expected to have unusual properties. 

A related set of issues stems-from a gedankenexper- 
iraent proposed by AndersonJ13 in which two previ- 
ously isolated Bose condensates are brought into particle- 
exchanging contact. Before contact is initiated, each 
isolated subsystem has its own macroscopically occu- 
pied state. The two isolated Bose subsystems consid- 
ered together therefore provide a simple example of a 
system with a "fragmented" condensate. What happens 
when the condensates are allowed to exchange particles 
- will the phase rdifierence be random from experiment 
to experiment JI9t2MiJ or be reproducibly zerotlrta? 

Here we show that when fragmented condensates over- 
lap, they are inherently unstable to the formation of a 
single, conventional condensate of well-defined relative 
phase. The essential idea is that even weak Josephson 
coupling between the fragments rapidly generates phase 
coherence, modifying the density matrix deterministi- 
cally to give a unique macroscopically occupied single 
particle state. The proof proceeds by first assuming the 
existence of a Bose state with multiple condensates, and 
then explicitly constructing a lower energy state. The 



argument is very general, and applies to strongly as well 
as weakly interacting systems, and for either sign of the 
scattering length. As a corollary, we show that the rela- 
tive phase that emerges when two previously uncoupled 
condensates are placed in contact is reproducible, and is 
determined by the physics of coherent particle exchange 
between them. 

Multiple condensates. Consider a fluid of N identi- 
cal Bose particles. Let |0) be a many-body state that is 
presumed to be both (a) the ground state of a many-body 
Hamiltonian H and (b) a "fragmented condensate," i.e., 
a state whose reduced one-body density matrix 



P(r, 



(^(r')^(r)) 



(1) 



has two or more macroscopic eigenvalues Ni, with cor- 
responding eigenfunctions 4>i{r). The Ni have the phys- 
ical interpretation of the mean occupancy of the single- 
particle state 4»i in the many-body state |0). In general 
the macroscopic Ni are not integers, and they need not 
add up to the total particle number N. For simplicity, 
we will assume that ony two such macroscopic eigenval- 
ues exist, but our arguments are easily generalized to any 
finite number a 

We define the condensate annihilation operators 



h= / dr&(r)V(r) 



(2) 



where is the particle annihilation operator at posi- 
tion r. Then in the fragmented state, 



[polij = (0\b]b i \0) = N i 5 ij , 



(3) 



which follows from the assumption that 4>i are eigenstates 
of p. For i 7^ j the density matrix element vanishes, indi- 
cating the absence of coherence between the condensates. 

A family of states. We now construct a family of N- 
particle states labeled by the integers q — 0, ±1, ±2, ... in 
which q particles are transferred between the two con- 
densates. For positive q, the transfer is from 1 to 2; for 
negative q, the transfer is from 2 to 1; the q = state 
is simply the original fragmented state |0). These states 
can be formed by repeated applications of the operator 



(4) 



for positive q, and its Hermitian conjugate B for negative 
q. The factors in the denominator make the states 



|g> = (i?t)«|0) q>0 
|ff) = (B)l«l|0) q<0 



(5) 
(6) 



normalized for q <C Ni , which will be the case of interest 
below. The operator Q can then be defined as 9k) (<? I- 
Our strategy will be to construct a linear combination of 
the states \q) that is lower in energy than the putative 
ground state |0). 
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An effective tight-binding model. We consider a 
many-body Hamiltonian TL given by the sum of a one- 
body part h Q - kinetic energy plus a single-particle po- 
tential V(r) - and a two-body interaction U(r—r'), which 
for notational simplicity we will take to be a (5-function 
pseudopotential of strength U . (Our results are easily 
extended to more general interactions.) Within the sub- 
space of states {!<?)}, H can be written 

H = -[tiB+ + t\B] - [t 2 (fit)2 + t * B 2} + 1 Q2 ; ( 7 ) 

which has the form of a one-dimensional tight-binding 
model with nearest and next-nearest neighbor hopping 
matrix elements ti and i 2 t and a harmonic potential of 
spring constant k. 

The generally complex hopping matrix elements are 
given by 



t x = -VN^N 2 J dr </)*(r)Wi(r), (8) 
t 2 = ~NiN 2 J dr^(r)^(r)0i(r)^i(r), (9) 

where h e g = h +U^2 i iVj|</>j(r)| 2 is an effective one-body 
Hamiltonian that includes mean field interactions. 

If the condensates </>i(r) extend over the same d- 
dimensional volume R d , where R is a characteristic 
length scale, then the matrix elements of eq. (|^) scale as 
t\ ~ Nt\ and t 2 ~ Ne 2 , where the are functions of the 
mean density N/R d with the units of energy per particle. 
Typically, they will be p£the same order of magnitude as 
the chemical potentialrj We consider below the special 
case in which the matrix elements t\ and t 2 vanish for 
symmetry reasons. 

The hopping terms of the Hamiltonian (fy give rise to 
a tight-binding band with dispersion 



E(k) 



[he 



-i2k 



+ t* 2 e l2k ] 



(10) 



which can be expanded about its minimum to yield 

TN T'N , N9 , x 

E(k) = -— + —-(k - x) 2 + ... (11) 

where the real numbers T and T' depend on the e il and 
X can be found by minimizing eq. ((l^). The effective 
mass for this band is given by ft 2 /M e s — T'N . 

The "spring constant" k in eq. (^) can be found quite 
generally from the assumption that the condensates are 
compressible, which is expected on general grounds for a 
Bose fluid. E£l Then the energy of a state with qi particles 
added to condensate i (by repeated application of b\) 
will be a smooth function that can be Taylor expanded 
for small qi <C Ni, so that 



(q\H\q) = E Q + qifn - fi 2 ) + -q 2 



(12) 



where /i, = dE/dNi is the chemical potential of the i-th 
fragment and 



K 
N 



r d 2 E 



d 2 E 



L dN{ 8N 1 dN 2 



d 2 E 
dN 2 



(13) 



For condensates <j>i of volume R d , if is a number of order 
U{N/R d ) that is again comparable to the chemical po- 
tential, and hence T and X". Since by assumption q = 
is the ground state the chemical potentials must be equal 
and the linear term in q vanishes. For stability, k must 
be positive. 

The tight-binding model ([?]) for two coupled conden- 
sates is related by a canonically transformation to the 
familiar pendulum description of a Josephson junction,Ej 
which treats the relative phase as a periodic "coordi- 
nate" and the number fluctuation as its corresponding 
"momentum." The present analysis instead takes the 
number fluctuation as a discrete coordinate, with rela- 
tive phase corresponding to the conjugate (crystal) mo- 
mentum. The two approaches are completely equivalent, 
and provide complementary insights to the behavior of 
coupled condensates. 

A lower energy state. Within a continuum approx- 
imation, the tight-binding Hamiltonian (Q) describes a 
particle of mass M e g in a harmonic potential of spring 
constant k. The frequency of this oscillator is given by 
M c sw 2 g = K, so that huj c ff — \]T'K. The ground state 
has the form 



\G) 



v ^ e 



<\q), 



with a characteristic spread in q 

a = (VMcff^eff) 1 / 2 - (T'/KY^N 1 / 2 , 



(14) 



(15) 



Since a ~ iV 1 / 2 , we are justified a posteriori in both 
the continuum approximation and our decision to neglect 
terms of order q 3 (and higher) in E(q). 
The energy of \G) relative to |0) is 



Eg — E 



TN (T'K) 1 / 2 
~~ + 2 : 



(16) 



where the first term on the right is the delocalization 
energy of the bottom of the tight-binding band, and the 
second is the zero-point energy of a harmonic oscillator 
with energy spacing ttu> c s- Evidently the admixture of 
correlated condensate fluctuations \G) is lower in energy 
than the fragmented state \q = 0) if N >-. T'K/(T 2 ), 
where the right hand side is of order unityEj Thus even 
for extremely weak off-diagonal couplings, a fragmented 
condensate will not be the ground state for a macroscopic 
number of particles. 

Density matrix revisited. What is the nature of 
the true ground state \G)1 It is easy to show that \G) 
is, in fact, a conventional Bose condensate whose density 
matrix has a unique macroscopic eigenvalue. Since q <ti 
Ni, the diagonal matrix elements of p are unchanged 



[p G ]u = {G\b\bi\G) = Ni. 



(17) 
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The off-diagonal matrix elements, however, are now 
macroscopic: 

[p G ]i2 = (G\bth\G) = g e-*Xy/NxN 2 (18) 

where 



where A is a number of order unity. 

Rediagonalizing the density matrix, we find a single 
macroscopically occupied eigenstate 

fa(r)= Cl fa(r)+e^c 2 fa(r), (20) 

where = y/Ni/[N\ + A^]. The occupation of this state 
is N g = Ni + N 2 — B, where B is proportional to A, and 
is also a number of order unity. The orthogonal combi- 
nation 4>-{r) = c 2 fa(r) — e lx ci4> 2 (r) has eigenvalue B of 
order unity. We conclude that fragmented condensates 
are unstable towards the formation of conventional Bosc 
condensate with a unique macroscopically occupied state. 

Spontaneous symmetry breaking. As noted 
above, the matrix elements may vanish for symme- 
try reasons. As a specific example, consider a rota- 
tionally invariant system with two fragments fa and fa 
that have angular momentum projections mi, m%, re- 
spectively. Then the ti will be zero. This situation 
arises in the context of rotating Bose gases trapped in 
an azimuthally symmetric potential.!!! Under these cir- 
cumstances, we must ask if the fragmented condensate 
remains stable in the presence of a small symmetry- 
breaking perturbation, which introduces a nonzero t%. 

We have seen, however, that even for T ~ \/N 2 a frag- 
mented condensate will be unstable to the formation of 
a conventional condensate ( p0| ) that is a superposition of 
the fragments, with a relative phase x that is determined 
by the details of the perturbing potential through eq. (||). 
The resulting state is not an eigenstate of angular mo- 
mentum, and therefore has an asymmetric density and 
current distribution.Ej The fact that an order 1/N per- 
turbation can reduce the symmetry of the ground state 
indicates that fragmented condensates will spontaneously 
break whatever symmetry (rotation, in the case of ref. |^) 
permitted fragmentation in the first place. 

This result can be understood in another way by con- 
sidering the susceptibility of a fragmented condensate to 
a perturbation which couples to the inter-fragment cur- 
rent density J = i[t x B^ -t\B\. Since ((B 1 ) 2 ) = (B 2 ) = 
in the q = state, while B^B = BB^ = 1 as an opera- 
tor identity, the mean-square current fluctuation ( J 2 ) ~ 
|<i | 2 . A weak one-body potential that allows scattering 
between fa and fa, however, will typically introduce a 
t\ r~j N. Fragmented condensates are thus highly sus- 
ceptible to perturbations that permit coherent particle 
transfer between the fragments. 

Gauge covariance. The phase factor e %x which en- 
ters into the superposition ( |20| ) is not arbitrary, but is 



determined by the phase of the off-diagonal density ma- 
trix element (|l8|), which in turn is given by the phase 
of the Hamiltonian matrix elements t\ and t 2 that scat- 
ter particles between the two condensates. We can easily 
confirm that this result is properly gauge covariant.Q 

The single-particle states fa are eigenstates of the 
(fragmented) one-body density matrix po, and are only 
defined to within arbitrary phase factors. All physical ob- 
servables must therefore be unchanged if each fa is mul- 
tiplied by e lai . It is easy to follow these phases through 
our calculations. The hopping matrix elements trans- 
form as ti -> t 1 e l{ai - a ^ and t 2 -> i 2 e 2i(Ql_Q2) , so that 
X — > X + a i ~ a 2- At the end of the day, the unique 
condensate ( pp| ) becomes 

fa(r) -» ci[e iQl 0i(r)] + e^ x+ ^- a ^[c 2 e la - fa(v)} 

= ^fa{v). (21) 

The arbitrary phases a; are only reflected in the overall 
phase of far (which itself is only defined up to an over- 
all phase). The relative phase of </> g (r) between any two 
points ri and r 2 are invariants. The superfluid veloc- 
ity, which is proportional to the gradient of the phase, is 
also invariant. Thus the phase of the off-diagonal one- 
body density matrix element (proportional to (B)) is a 
gauge-covariant representation of the relative phase be- 
tween two condensates. 

Relative phase of previously isolated conden- 
sates. We can use a similar approach to analyse Ander- 
son's gedankenexperimenMl in which two independently 
created condensates 4>l and 4>r are placed in contact. 
Before contact is established, the many-body state of the 
combined system is a simple product of states - i.e., a 
"fragmented condensate" - that we can denote \q = 0). 
Since the off-diagonal element (q = 0|6^6^|g = 0) of the 
density matrix is zero, the relative phase between the two 
condensates is the phase of the complex number zero, 
which is ill-defined. 

For gedtmfcenpurposes, let us assume that the two in- 
dependent condensates are identical, with equal chemi- 
cal potentials. (Our discussion is easily generalized.) As 
shown in the previous section, we may without loss of 
generality choose both fa^ and fan to be real. The time 
evolution following the initiation of contact is goiierned 
by a Hamiltonian of the form (Q), plus damping,Ej with 
i 2 << t\. By time-reversal symmetry, the are real, 
and x — or 7r; for definiteness we assume conventional 
tunneling contact, with \ = 0. In contrast to the case 
of spatially overlapping fragments, however, the Joseph- 
son coupling t c g between two reservoirs in contact will 
be proportional to their contact area, and-exponentially 
small in the energy barrier between them.tS 

Except for the discreteness of q, this is a familiar prob- 
lem. At long times, the system ends up in a Gaussian 
ground state, which we have seen is simply a conven- 
tional condensate with uniform phase. Relative to this 
"vacuum," the initial state q — is "squeezed," with 
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a width 6q ~ 1 in position and 8p ~ it in the conju- 
gate momentum. The establishment of phase coherence 
between two previously isolated condensates is therefore 
equivalent to. the decay of a "squeezed vacuum" in quan- 
tum opticsEj For our purposes, the important feature is 
that due to the ±q symmetry of the initial state and the 
Hamiltonian, (B^) = (B) for all time. The off-diagonal 
density matrix element is then always real, so as soon as 
the relative phase becomes denned,_Lt is zero. 

The experiment of Andrews et alE3 in which two nom- 
inally identical, initially isolated condensates are allowed 
to expand and overlap can also be cast in this form. If 
the condensates are initially converging symmetrically, 
then 4>r(z) — 4>i,(—z)*, and we again find that the 
transfer matrix element that scatters particle between 
the condensates is real. This, result is consistent with 
Naraschewski et aL's analysiscj of the interference fringes 
seen in the MIT experiment, which requires \ — 0. 

Random vs. reproducible phases. Javanainen and 
YocO and Castin and Dalibardlla have discussed proto- 
cols in which a series of measurements is performed on 
the states of individual particles removed from two iso- 
lated condensates. Their elegant analyses show that the 
interference patterns that emerge are consistent with a 
well-defined phase difference between the two conden- 
sates. This phase difference, however, is arbitrary, and 
varies if the entire series of measurements are repeated. 
How can our calculation be reconciled with these results? 

In the schemes of refs. [l3] and [l4], each measurement 
of the series removes a particle from the two-condensate 
system in a defined superposition of the original conden- 
sates (pi. The remaining particles are therefore left in 
a specific entangled many-body state that depends on 
the results of the measurements. In other words, a se- 
ries of observations of the removed particles generates 
a specific perturbation 5TL(t) acting on the remaining 
particles. This perturbation does not conserve particle 
number, and in effect applies,a time-dependent, gauge- 
symmetry-breaking "r/-field"E3 to the condensates. 

As we have seen, the susceptibility of the state q = 
to such a perturbation is divergent with N. Thus it only 
takes a small perturbation (i.e., a few measurements) to 
sculpt the state of the remaining particles into a single 
condensate, whose relative phase depends deterministi- 
cally on the STC(t), and hence on the specific sequence 
of observations performed on the removed particles. We 
have shown above that when the two condensates are 
placed in particle-exchanging contact their coupling also 
introduces a perturbation, which reproducibly enforces 
a relative phase of zero. An interesting question is the 
manner in which this effect competes with the measure- 
ment scenarios of refs. |l3| and |i~4| . 
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the unique ground state can be determined. This many- 
body ground state is a conventional Bose condensate with 
a unique condensate wavefunction. The crucial requirement 
is that the occupation numbers of the condensates must be 
macroscopic, hence the restriction to a finite number of con- 
densates. It is an interesting question to ask whether the 
number of condensates N con can scale in some nontrivial 
way with total particle number N to preserve fragmenta- 
tion. 

We emphasize that the behavior of multiple spatially over- 
lapping condensates (as proposed by Ref. ^, and considered 
in the first part of the present work) can be different from 
that of two spatially separated condensates (whose behav- 
ior when brought into contact is discussed in the second 
part of this paper). For separated condensates the hopping 
matrix elements ti will be exponentially small in the height 
of the intervening energy barrier. Then TN, while formally 
macroscopic (i.e., proportional to iV), can also be negligi- 
bly small. For overlapping condensates this behavior is not 
expected. 

The assumption that the condensate fragments are com- 
pressible fluids excludes the case of Bose insulators in 
particle-exchanging contact, for which the energy (q\H\q) 
would depend on \q\. 

We emphasize again that for the case of two condensates 
separated by a barrier, T and T' will be exponentially small 
in the barrier height. Thus according to the condition de- 
rived here, a "fragmented" state could then be stable even 
for "macroscopic" N. This is simply the statement that two 
sufficiently far separated reservoirs of bosons can condense 
independently. This behavior is distinct from the fragmen- 
tation of spatially overlapping condensates as discussed in 
ref. §. 



